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Abstract 

Using the renormalization group method, new type of fluctuation-driven first 
order phase transitions and critical phenomena are predicted for certain classes of 
ferromagnetic superconductors and superfluids with unconventional (spin-triplet) 
Cooper pairing. The problem for the quantum phase transitions at extremely low 
and zero temperatures is also discussed. The results can be applied to a wide class 
of ferromagnetic superconductive and superfluid systems, in particular, to itinerant 
ferromagnets as UGe2 and URhGe. 
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1. Introduction 

In this paper an entirely new critical behavior in unconventional ferromagnetic superconductors 
and superfluids is established and described. This phenomenon corresponds to an isotropic 
ferromagnetic order in real systems but does not belong to any known universality class pQ 
and, hence, it might be of considerable experimental and theoretical interest. Due to crystal 
and magnetic anisotropy a new type of fluctuation-driven first order phase transitions occur, 
as shown in the present investigation. These novel fluctuation effects can be observed near 
finite and zero temperature ( "quantum" ) phase transitions [U [2] in a wide class of ferromagnetic 
systems with unconventional (spin-triplet) superconductivity or superfluidity. 

The present investigation has been performed on the concrete example of intermetallic com- 
pounds UGe2 and URhGe, where the remarkable phenomenon of coexistence of itinerant fer- 
romagnetism and unconventional spin-triplet superconductivity [3] has been observed [4]. For 
example, in UGe2, the coexistence phase occurs [4] at temperatures < T < 1 K and pres- 
sures 1 < P < Pq ~ 1.7 GPa. A fragment of (P,T) phase diagrams of itinerant ferromagnetic 
compounds (5j is sketched in Fig. 1, where the lines Tp{P) and T C {P) of the paramagnetic (P) - 
to-ferromagnetic(F) and ferromagnetic-to-coexistence phase(C) transitions are very close to each 
other and intersect at very low temperature or terminate at the absolute zero (PoiO). At low 
temperature, where the phase transition lines are close enough to each other, the interaction 
between the real magnetization vector M{r) = {Mj(r);j = 1, ...,m} and the complex order pa- 
rameter vector of the spin-triplet Cooper pairing [3], ip(r) = {ijj a (r) = {ip' a -\-iipa)) a = I, ....n/2} 
(n = 6) cannot be neglected pQ and, as shown here, this interaction produces new fluctuation 
phenomena. 
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Fi gure 1: (P,T) diagram with a zero-temperature multicritical point (Po,0). Para- (P), ferromagnetic 
(F), and coexistence (C) phases, separated by the lines Tf(P) and T C (P) of P-F and F-C phase transitions, 
respectively. 

Both thermal fluctuations at finite temperatures (T > 0) and quantum fluctuations (correla- 
tions) near the P-driven quantum phase transition at T = should be considered but at a 
first stage the quantum effects [2] can be neglected as irrelevant to finite temperature phase 
transitions (Tp ~ T c > 0) . The present treatment of a recently derived free energy functional [3] 
by the standard Wilson-Fisher renormalization group (RG) [I] shows that unconventional ferro- 
magnetic superconductors with an isotropic magnetic order (m = 3) exhibit a quite particular 
multi-critical behavior for any T > 0, whereas the magnetic anisotropy (m = 1,2) generates 
fluctuation-driven first order transitions pp. Thus the phase transition properties of spin-triplet 
ferromagnetic superconductors are completely different from those predicted by mean field theo- 
ries [5, 6j. The results can be used in the interpretation of experimental data for phase transitions 
in itinerant ferromagnetic compounds [7j- 

The study presents for the first time an example of complex quantum criticality characterized 
by a double-rate quantum critical dynamics. In the quantum limit (T — > 0) the fields M and 
ip have different dynamical exponents, zm and z^, and this leads to two different upper critical 
dimensions: dff = 6 — zm and d¥ = 6 — z^. The complete consideration of the quantum 
fluctuations of both fields M and ip requires a new RG approach in which one should either 
consider the difference (zm — z^) as an auxiliary small parameter or create a completely new 
theoretical paradigm of description. The considered problem is quite general and presents a 
challenge to the theory of quantum phase transitions [2j. The results can be applied to any 
natural system within the same class of symmetry although this report is based on the example 
of itinerant ferromagnetic compounds. 

2. Renormalization-group investigation 

The relevant part of the fluctuation Hamiltonian of unconventional ferromagnetic superconduc- 
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tors [B] can be written in the form 



H = £ (r + A; 2 ) ^(k)\ 2 + i ( f + k 2 ) \M{k)\A + 
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M (fci) . [V (fe 2 ) x V* (fci + fca)] 



(1) 



where V ~ L d is the volume of the d— dimensional system, the length unit is chosen so that 
the wave vector k is confined below unity (0 < k = \k\ < 1), g > is a coupling constant, 
describing the effect the scalar product of M and the vector product (ip x ip*) for symmetry 
indices m = (n/2) = 3, and the parameters t ~ (T — Tf) and r ~ (T — T s ) are expressed by the 
critical temperatures of the generic (g = 0) ferromagnetic and superconducting transitions. As 
mean field studies indicate [HE], T S (P) is much lower than T C (T) and Tp(P) / Tf(P). 

The fourth order terms (M 4 , |?/>| 4 , M 2 ^ 2 ) in the total free energy (effective Hamiltonian) [5j[6] 
have not been included in Eq. (1) as they are irrelevant to the present investigation. The simple 
dimensional analysis shows that the g— term in Eq. (1) corresponds to a scaling factor & 3_d / 2 and, 
hence, becomes relevant below the upper borderline dimension djj = 6, while fourth order terms 
are scaled by a factor b 4 ~ d as in the usual 4 — theory and are relevant below d < 4 (b > 1 is a 
scaling number) pQ. Therefore we should perform the RG investigation in spatial dimensions d = 
6 — e where the g-term in Eq. (1) describes the only relevant fluctuation interaction. Moreover, 
the total fluctuation Hamiltonian [5l[6] contains off-diagonal terms of the form kikji/j a tpt; i ^ j 
and/or a ^ (3. Using a convenient loop expansion these terms can be completely integrated out 
from the partition function to show that they modify the parameters (r, t, g) of the theory but 
they do not affect the structure of the model (1). So, such terms change auxiliary quantities, 
for example, the coordinates of the RG fixed points (FPs) but they do not affect the main RG 
results for the stability of the FPs and the values of the critical exponents. Here we ignore these 
off-diagonal terms. 

One may consider several cases: (i) uniaxial magnetic symmetry, M = (0, 0, M3), (ii) tetragonal 
crystal symmetry when ip = (ipi, ip2i 0)> (hi) XY magnetic order (Mi, M2, 0), and (iv) the general 
case of cubic crystal symmetry and isotropic magnetic order (m = 3) when all components of the 
three dimensional vectors M and tp may have nonzero equilibrium and fluctuation components. 
The latter case is of major interest to real systems where fluctuations of all components of 
the fields are possible despite the presence of spatial crystal and magnetic anisotropy that 
nullifies some of the equilibrium field components. In one-loop approximation, the RG analysis 
reveals different pictures for anisotropic (i)-(iii) and isotropic (iv) systems. As usual, a Gaussian 
("trivial") FP (g* = 0) exists for all d > and, as usual pQ this FP is stable for d > 6 where the 
fluctuations are irrelevant. In the reminder of this paper the attention will be focussed on spatial 
dimensions d < 6, where the critical behavior is usually governed by nontrivial FPs (<?* 7^ 0). In 
the cases (i)-(iii) only negative ( "unphysical" |9j) FP values of g 2 have been obtained for d < 6. 
For example, in the case (i) the RG relation for g takes the form 



where g' is the renormalized value of g, r/ = (K c [-\/8)g 2 is the anomalous dimension (Fisher's 
exponent) [I] of the field M 3 ; K d = 2 1 - d -K- d / 2 /Y(d/2). Using Eq. (2) one obtains the FP 
coordinate (g 2 )* = — 967r 3 e. For d < 6 this FP is unphysical and does not describe any critical 
behavior. For d > 6 the same FP is physical but unstable towards the parameter g as one may 
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see from the positive value y g = — lle/2 > of the respective stability exponent y g defined 
by &9 1 = t> V9 5g. Therefore, a change of the order of the phase transition from second order in 
mean- field ( "fluctuation free" ) approximation to a fluctuation-driven first order transition when 
the fluctuation ^-interaction is taken into account takes place. This conclusion is supported by 
general concepts of RG theory [1] and by the particular property of these systems to exhibit 
first order phase transitions [6] in mean field approximation for broad variations of T and P. 

In the case (iv) of isotropic systems the RG equation for g is degenerate and the e-expansion 
breaks down. A similar situation is known from the theory of disordered systems [9] but here 
the physical mechanism and details of description are different. Namely for this degeneration 
one should consider the RG equations up to the two-loop order. The derivation of the two-loop 
terms in the RG equations is quite nontrivial because of the special symmetry properties of the 
interaction g-term in Eq. (1). For example, some diagrams with opposite arrows of internal 
lines, as the couple shown in Fig. (2), have opposite signs and compensate each other. The 
terms bringing contributions to the (/-vertex are shown diagrammatically in Fig. 3. The RG 
analysis is carried out by a completely new e 1//4 -expansion for the FP values and e 1//2 -expansion 
for the critical exponents; again e = (6 — d). The RG equations are quite lengthy and here only 
the equation for g is discussed. It has the form 



g' = b( e ~ 2 ^-VM)/2 g j-j + Ag 2 + 3 ^ 2B + c y 



(3) 
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ijm and r/if, are the anomalous dimensions of the fields M and ift, respectively. The one-loop 
approximation gives correct results to order e 1 / 2 and the two-loop approximation brings such 
results up to order e. In Eq. (4), r and t are small expansion quantities with equal FP values 
t* = r* = Kdg 2 . Using the condition for invariance of the two fc 2 -terms in Eq. (1) one obtains 
f]M = Vip = V, where 
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Eq. (3) yields a new FP 
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which corresponds to the critical exponent r\ = 2(2e/13) 1 / 2 — 2e/3 (for d = 3, 
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The eigenvalue problem for the RG stability matrix M. = [(dfMi/dfj,j); (/ii, /X2, M3) 
be solved by the expansion of the matrix elements up to order e 3//2 . When the eigenvalues 
Xj = Aj(b)b y J of M are calculated dangerous large terms of type b 2 and 6 2 (ln6), (b 1) [8] in 
the off-diagonal elements of the matrix M ensure the compensation of redundant large terms 
of the same type in the diagonal elements Ma. This compensation is crucial for the validity of 
scaling for this type of critical behavior. Such a problem does not appear in standard cases of 
RG analysis [U |8] . As in the usual c/> 4 -theory [8] the amplitudes Aj depend on the scaling factor 
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Figure 2: A sum of <7 5 -diagrams equal to zero. The thick and thin lines correspond to correlation 
functions (\ip a \ 2 ) and (|Mj| 2 ), respectively; vertices (•) represent g-term in Eq. (1). 

b: Ai = A 2 = 1 + (27/13)6 2 e, A 3 = 1 - (81/52)e(ln6) 2 . The critical exponents y t = y x , y r = y 2 
and y g = y% are 6-invariant: 



2e 197 

Vr = 2 + 10a/ 1 e. 

y V 13 39 ' 



(9) 



Vt = Vr ~ 18(2e/13) 1//2 , and y g = — e > for d < 6. The correlation length critical exponents 
u$ = l/y r and vm = 1/yt corresponding to the fields ip and M are 



1 5 2e 103 

2 ~ 2 V 13 + 156 6 ' 



(10) 
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These exponents describe a quite particular multi-critical behavior which differs from the nu- 
merous examples known so far. For d = 3, = 0.78 which is somewhat above the usual value 
v ~ 0.6 ~r~ 0.7 near a standard phase transition of second order [I], but i^m = 1-76 at the same 
dimension (d = 3) is unusually large. The fact that the Fisher's exponent rj is negative for 
d = 3 does not create troubles because such cases are known in complex systems, for example, 
in conventional superconductors [10]. Perhaps, a direct extrapolation of the results from the 
present e-series is not completely reliable because of the fact that the series has been derived 
under the assumptions of e -C 1 and under the conditions e l l 2 b < 1, e 1 / 2 (lnfe) -C 1 provided 
b > 1. These conditions are stronger than those in the usual </) 4 -theory [U[8]. Using the known 
relation PQ 7 = (2 — r])u, the susceptibility exponents for d = 3 take the values 7^, = 2.06 and 
7m = 4.65. These values exceed even those corresponding to the Hartree approximation [l] 
(7 = 2v = 2 for d = 3) and can be easily distinguished in experiments. Note, that here we 
follow the interpretation of the asymptotic e-series in the way given by Lawrie et al. [9]. This 
point of view is quite comprehensive, in particular, for an avoiding artificial conclusions from 
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Figure 3: Diagrams for g' of third and fifth order in g. The arrows of the thick lines have been omitted. 

the RG analysis of complex systems with competing effects, such as the systems described by 
the Eq. (1). 

Notes about the quantum effects on the phase transitions. The critical behavior discussed so 
far may occur in a close vicinity of finite temperature multi-critical points (T c = Tj > 0) in 
systems possessing the symmetry of the model (1). In certain systems, as shown in Fig. 1, this 
multi-critical points may occur at T = 0. In the quantum limit (T — * 0), or, more generally, 
in the low-temperature limit [T <C \i\ \i = (t,r);ks = 1] the thermal wavelengths of the fields 
M and ip exceed the inter-particle interaction radius and the quantum correlations fluctuations 
become essential for the critical behavior OUT]. The quantum effects can be considered by RG 
analysis of a comprehensively generalized version of the model (1), namely, the action S of the 
referent quantum system. The generalized action is constructed with the help of the substitution 
(—Tt/T) — > S[M(q),ip(q)]. Now the description is given in terms of the (Bose) quantum fields 
M(q) and i/j(q) which depend on the (d + l)-dimensional vector q = (u>i,k); ui = 2-irlT is the 
Matsubara frequency (H = 1; I = 0, ±1, . . . ). The fc-sums in Eq. (1) should be substituted by 
respective g-sums and the inverse bare correlation functions (r + k 2 ) and (t + k 2 ) in Eq. (1) 
contain additional u;;— dependent terms, for example [21 [TT] 

(\i> a (q)\ 2 )~ l = \"l\ + k 2 + r. (12) 

The bare correlation function (\Mj(q)\) 2 contains a term of type \uJi\/kP , where 6 = 1 and 9 = 2 
for clean and dirty itinerant ferromagnets, respectively [TT]. The quantum dynamics of the field 
tp is described by the bare value z = 2 of the dynamical critical exponent z = zw, whereas the 
quantum dynamics of the magnetization corresponds to zm = 3 (for 9 = 1), or, to zm = 4 
(for 9 = 2). This means that the classical-to-quantum dimensional crossover at T — > is given 
by d — > (d + 2) and, hence, the system exhibits a simple mean field behavior for d > 4. Just 
below the upper quantum critical dimension d^ = 4 the relevant quantum effects at T = 
are represented by the field tp whereas the quantum (w/-) fluctuations of the magnetization are 
relevant for d < 3 (clean systems), or, for even for d < 2 (dirty limit) [11] . This picture is 
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confirmed by the analysis of singularities of the relevant perturbation integrals. Therefore the 
quantum fluctuations of the field ip have a dominating role below spatial dimensions d < 4, and 
for dimensions 3 < d < 4 (clean systems), or, for 2, d < 4 in case of dirty limit, they are the only 
quantum fluctuations in these systems. 

Taking into account the quantum fluctuations of the field ip and completely neglecting the in- 
dependence of the magnetization M , cq = (4 — (f)-analysis of the generalized action S has been 
performed within the one-loop approximation (order eg). In the classical limit (r/T -C 1) one 
re-derives the results already reported above together with an essentially new result, namely, 
the value of the dynamical exponent 2W> = 2 — (2e/13) 1//2 which describes the quantum dynamics 
of the field if). In the quantum limit (r/T ^> 1, T — * 0) the static phase transition properties 
are affected by the quantum fluctuations, in particular, in isotropic systems (n/2 = m = 3). 
For this case, the one-loop RG equations corresponding to T = are not degenerate and give 
definite results. The RG equation for g, 



yields two FPs: (a) a Gaussian FP (g* = 0), which is unstable for d < 4, and (b) a FP 
(g 2 )* = — 127r 3 eo which is unphysical [(g 2 )* < 0] for d < 4 and unstable for d > 4. Thus the 
new stable critical behavior corresponding to T > and d < 6 disappears in the quantum limit 
T — > 0. At the absolute zero and any dimension d > the P— driven phase transition (Fig. 1) 
is of first order. This can be explained as a mere result of the limit T — > 0. The only role of the 
quantum effects is the creation of the new unphysical FP (b). In fact, the referent classical system 
described by H from Eq. (1) also looses its stable FP (8) in the zero-temperature (classical) limit 
T — ► but does not generate any new FP because of the lack of g 3 -term in the equation for 
g'; see Eq. (13). At T = the classical system has a purely mean field behavior [2] which is 
characterized by a Gaussian FP (g* = 0) and is unstable towards T-perturbations for < d < 6. 
This is a usual classical zero temperature behavior where the quantum correlations are ignored. 
For the standard <p 4 — theory this picture holds for d < 4. One may suppose that the quantum 
fluctuations of the field ip are not enough to ensure a stable quantum multi-critical behavior at 
T c = Tp = and that the lack of such behavior in result of neglecting the quantum fluctuations 
of M. One may try to take into account these quantum fluctuations by the special approaches 
from the theory of disordered systems, where additional expansion parameters are used to ensure 
the marginality of the fluctuating modes at the same borderline dimension djj (see, e.g., Ref. [S]). 
It may be conjectured that the techniques known from the theory of disordered systems with 
extended impurities cannot be straightforwardly applied to the present problem and, perhaps, 
a completely new supposition should be introduced. 

3. Final remarks 

The present results may be of use in interpretations of recent experiments [7] in UGe2, where 
the magnetic order is uniaxial (Ising symmetry) and the experimental data, in accord with the 
present consideration, indicate that the C-P phase transition is of first order. Systems with 
isotropic magnetic order are needed for an experimental test of the new multi-critical behavior. 
Besides, the present investigation exhibits several new essential problems which are a challenge 
to the theory of quantum phase transitions. 
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